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ABSTRACT 



Aims. Results of the long-term (11 years, from 1996 to 2006) Ho- and HjS line variations of the active galactic nucleus of NGC 4151 
are presented. 

Methods. High quality spectra (S/N> 50 and R«8 A) of Ho- and HyS were investigated. We analyzed line profile variations during 
monitoring period. Comparing the line profiles of Hor and Hj8, we studied difi'erent details (bumps, absorbtion features) in the line 
profiles. The variations of the different Hor and H/3 line profile segments have been investigated. Also, we analyzed the Balmer 
decrement for whole line and for line segments. 

Results. We found that the line profiles were strongly changing during the monitoring period, showing blue and red asymmetries. 
This indicates a complex BLR geometry of NGC 4151 with, at least, three kinematically distinct regions: one that contributes to the 
blue line wing, one to the line core and one to the red line wing. Such variation can be caused by an accelerating outflow starting very 
close to the black hole, where the red part may come from the region closer to the black hole than the blue part, which is coming 
from the region having the highest outflow velocities. 

Conclusions. Taking into account the fact that the BLR of NGC 4151 has a complex geometry (probably affected by an outflow) and 
that a portion of the broad line emission seems to have not a pure photoionization origin, one can ask the question whether the study 
of the BLR by reverberation mapping may be valid in the case of this galaxy. 
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1. Introduction 

In spite of many papers devoted to the physical properties 
(physics and geometry, see e.g. Sulentic et al. 2000) of the broad 
line region (BLR) in active galactic nuclei (AGN), the true na- 
ture of the BLR is not well known. The broad emission lines, 
their shapes and intensities can give us many information about 
the BLR geometry and physics. In the first instance, the change 
in the line profiles and intensities could be used for investigating 
the BLR nature. It is often assumed that variations of line profiles 
on long time scales are caused by dynamic evolution of the BLR 
gas, and on short time scales by reverberation effects (Sergeev et 
al. 2001). In a number of papers it was shown that individual seg- 
ments in the line profiles change independently on both long and 
short time scales (e.g.. Wanders and Peterson 1996; Kollatschny 
and Dietrich 1997; Newman et al. 1997; Sergeev et al. 1999). 
Moreover, the broad line shapes may reveal some information 
about the kinematics and structure of the BLR (see Popovic et 
al. 2004). 
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One of the most famous and best studied Seyfert galaxies 
is NGC 4151 (see e.g. UMch 2000; Sergeev et al. 2001; Lyuty 
2005, Shapovalova et al. 2008 - Paper I, and reference therein). 
This galaxy, and its nucleus, has been studied extensively at all 
wavelengths. The reverberation investigation indicates a small 
BLR size in the center of NGC 4151 (see e.g. Peterson and Cota 
1988: 6 + 4 l.d.; Clavel et al. 1990: 4 + 3 l.d.; Maoz et al. 1991: 
9 + 2 l.d.; Bentz et al. 2006: 6.6;^'^^ l.d.). 

Spectra of NGC 4151 show a P Cygni Balmer and He I ab- 
sorption with an outflow velocity from -500 km s"' to -2000 
kms ', changing with the nuclear flux (Anderson and Kraft 
1969; Anderson 1974; Sergeev et al. 2001; Hutchings et al. 
2002). This material is moving outward along the line of sight, 
and may be located anywhere outside ~15 light-days (UMch 
and Horne 1996). An outflow is also seen in higher velocity 
emission-line clouds near the nucleus (Hutchings et al. 1999), 
such as multiple shifted absorption lines in C IV and other UV 
resonance lines (Weymann et al. 1997; Crenshaw et al. 2000), 
while warm absorbers are detected in X-ray data (e.g. Schurch 
and Warwick 2002). 

Some authors assumed that a variable absorption is respon- 
sible, at least partially, for the observed continuum variability of 
AGN (Coflin-Souffrin et al.l996; Boiler et al. 1997; Brandt et 



al. 1999; Abrassart and Czerny 2000; Risaliti, Elvius, Nicastro 
2002). Czerny et al. (2003) considered that most of variations 
are intrinsic to the source, though a variable absorption cannot 
be quite excluded. 

The nucleus of NGC 4151 emits also in the radio range. The 
radio image reveals a 0.2 pc two-sided base to the well-known 
arc-second radio jet (Ulvestad et al. 2005). The apparent speeds 
of the jet components relative to the radio AGN are less than 
0.050c and less than 0.028c at nuclear distances of 0.16 and 6.8 
pc, respectively. These are the lowest speed limit yet found in a 
Seyfert galaxy and indicates non-relativistic jet motions, possi- 
bly due to thermal plasma, on a scale only an order of magnitude 
larger than the BLR (Ulvestad et al. 2005). 

The observed evolution of the line profiles of the Balmer 
lines of NGC 4151 was studied by Sergeev et al.(2001)in 1988- 
1998 and was well modeled within the framework of the two- 
component model, where two variable components with fixed 
line profiles (double-peaked and single-peaked) were used. 

Although the AGN of NGC 4151 has been much observed 
and discussed, there are still several questions concerning the 
BLR kinematics (disk, jets or more complex BLR) and dimen- 
sions of the innermost region. On the other hand, as was men- 
tioned above, multi-wavelength observations suggest both the 
presence of an accretion disk (with a high inclination) and an 
outflow emission (absorption). Consequently, further investiga- 
tions of the NGC 4151 nucleus are needed in order to constraint 
the kinematics, dimensions and geometry of its BLR. 

This work is subsequent to Paper 1 with the aim of study- 
ing the variations of both the integrated profiles of the broad 
emission lines and segments along the line profiles, during the 
(1 1-year) period of monitoring of NGC 4151. 

The paper is organized as follows: in §2 observations and 
data reduction are presented. In §3 we study the averaged spec- 
tral line profiles (over years, months and Periods I-Ul, see Paper 
I) of Ha and H/3, the line asymmetries and their FWHM vari- 
ations, light curves of different line segments and the line- 
segment to line-segment flux and continuum-line-segment re- 
lations. In §4 we analyze the Balmer decrements. The results are 
discussed in §6 and in §7 we outline our conclusions. 

2. Observations and data reduction 

Optical spectra of NGC 4151 were taken with the 6-m and 
1-m telescopes of SAO, Russia (1996-2006), with the 2.1- 
m telescope of the Guillermo Haro Astrophysical Observatory 
(GHAO) at Cananea, Sonora, Mexico (1998-2006), and with 
the 2. 1-m telescope of the Observatorio Astronomico Nacinal 
at San Pedro Martir (OAN-SMP), Baja California, Mexico 
(2005-2006). They were obtained with a long-slit spectrograph 
equipped with CCDs. The typical wavelength range was 4000 - 
7500 A, the spectral resolution was R=5-15 A, and the S/N ra- 
tio was > 50 in the continuum near Ho' and Hfi. In total 180 blue 
and 137 red spectra were taken during 220 nights. 

Spectrophotometric standard stars were observed every 
night. The spectrophotometric data reduction was carried out ei- 
ther with the software developed at the SAO RAS by Vlasyuk 
(1993), or with IRAF for the spectra observed in Mexico. The 
image reduction process included bias subtraction, flat-field cor- 
rections, cosmic ray removal, 2D wavelength linearization, sky 
spectrum subtraction, addition of the spectra for every night, 
and relative flux calibration based on standard star observations. 
Spectra were scaled to the constant flux F([0 5007). More 
details about observations and data reduction are given in Paper 
1 and will not be repeated. 



The observed fluxes of the emission lines were corrected for 
the position angle (PA), seeing and aperture effects (see Paper 1). 
The mean error (uncertainty) in our integral flux determinations 
for Ha and li/3 and for the continuum is <3%. In order to study 
the broad components of emission lines showing the main BLR 
characteristics, we removed from the spectra the narrow com- 
ponents of these lines and the forbidden lines. To this purpose, 
we construct spectral templates using the blue and red spectra in 
the minimum activity state (May 12, 2005). Both the broad and 
narrow components of H/3 and Ha, were fitted with Gaussians 
(see Fig. [1] available electronically only). The template spec- 
trum contains the following lines: for H/3 the narrow component 
of H/3 and [Oiii] AA 4959, 5007; for Ha the narrow component 
of Ha, [Nil] ^^6548, 6584, [Oi]/l^6300, 6364, [Sii]AA6in, 
6731. Then, we scaled the blue and red spectra according to our 
scaling scheme (see Appendix in Shapovalova et al. 2004), us- 
ing the template spectrum as a reference. The template spectrum 
and any observed spectrum are thus matched in wavelength, re- 
duced to the same resolution, and then the template spectrum is 
subtracted from the observed one. More details can be found in 
Paper 1 and Shapovalova et al. (2004). 

3. Line profile variations 

To investigate the broad line profile variations, we use the most 
intense broad lines in the observed spectral range, i.e. Ha and 
H/3, only from spectra with the spectral resolution of ~ 8 A. In 
Paper 1 we defined 3 characteristic time periods (1: 1996-1999, 
II: 2000-2001, III: 2002-2006) during which the line profiles of 
these lines were similar Average values and rms profiles of both 
lines were obtained for these periods. 

Here we recall some of the most important results of Paper 
I. In the first period (I, 1996-1999, JD=2450094.5-2451515.6), 
when the lines were the most intense, a highly variable blue com- 
ponent was observed, which showed two peaks or shoulders at 
~ -4000 km s ' and ~ -2000 km s ' in the rms Ha profiles 
and, to a less degree, in H/fl In the second period (II, 2000- 
2001, JD=245 1552.6-2452238.0) the broad fines were much 
fainter; the feature at ~ -4000 km s"' disappeared from the 
blue part of the rms profiles of both lines; only the shoulder at 

2000 km s"' was present. A faint shoulder at ~ 3500 km s"' 

was present in the red part of rms line profiles (see Fig. 6 in 
Paper I). In the third period (III, 2002-2006, JD=2452299.4- 
2453846.4) a red feature (bump, shoulder) at ~ 2500 km s"' 
was clearly seen in the red part of both the mean and the rms 
line profiles (see Fig. 7 in Paper I). In this paper we study the 
variations of the broad line profiles in more details. 

3.1. Month- and year-averaged profiles of the broad Ha and 
H/3 lines 

A rapid inspection of spectra shows that the broad line profiles 
vary negligibly within a one-month interval. On the other hand, 
in this time-interval, a slight variation in the broad line flux is 
noticed (usually around ~ 5 - 10%, except in some cases up to 
30%). Therefore we constructed the month-averaged line pro- 
files (see Fig.|2] available electronically only) of Ha and Hp. 

Moreover, the broad Ha and HyS line profiles were not chang- 
ing within a one-year period (or even during several years), while 
the broad fine fluxes sometimes varied by factors ~ 2-2.5 even 

' Here and after in the text the radial velocities are given with respect 
to the corresponding narrow components of Ha or H/3, i.e. it is accepted 
that V, = for the narrow components of Ha and H/3 
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Fig. 3. The year-averaged profiles (solid line) and their rms (dashed Une) of the Ha and H/? broad emission Unes in 1996-2006. The 
abscissae (OX) shows the radial velocities relative to the narrow component of the Ha or H/3 line. The ordinate (OY) shows the flux 
in units of 10"'"* erg cm"- s"' A"' . 



during one year. The smallest flux variations (factors ~ 1.1-1.3) 
were observed in 1996-1998 (during the line flux maximum). 
The largest line flux variations were observed in 2000-2001 and 
2005 (factors ~ 1.7 - 2.5), during the minimum of activity. 

As it was mentioned in Paper I, specific line profiles are ob- 
served during the three periods, but a more detailed inspection 
of the line profiles shows that sUght changes can also be seen be- 



tween the year-averaged profiles (see Fig. |3]l. Note here that in 
the central part of the Ho- profiles one can often see considerable 
residuals (e.g., in 2001 in the form of peaks) due to a bad subtrac- 
tion of the bright narrow components of Ha and [Nii] AA654E, 
6584, (at V, ~ -680 km s"' and ~ 960 km s"'). Therefore, we 
cannot conclude on the presence of some absorption in the cen- 
tral part of Ha. But, in the li/3 profiles (in the central part, at 



NGC4151 Hbeta NGC4151 Halphti 




-GOOD -4000 -2000 2DDD AuDD 600D -6000-4D0D-2000 2000 4000 6000 

Velocity [km s"'] Vsloclty [km s"'] 



Fig. 4. Some examples of month-averaged profiles of the Ha and 
Hj6 broad emission lines from 2000 to 2006. The abscissae shows 
the radial velocities relative to the narrow component of Ha or 
Hj6. The vertical dashed lines correspond to radial velocities: - 
2600 km s ' , km s ' and 3000 km s ' . The profiles are shifted 
vertically by a constant value. 



Vr ~ -430 kms"'; ~ -370 kms"') an absorption, especially 
strong from June 1999 to the end of 2000 (see Fig. |2]i, was de- 
tected. Note here that Hutchings et al. (2002) also found absorb- 
tion at yr(-2000 - 0) km s"' in the H/3 line. 

Now let us point some noticeable features in month-averaged 
line profiles: 

1) In 1996-2001 a blue peak (bump) in H/3 and a shoulder in 

Ha were clearly seen at 2000 km s ' (Fig. |2]i. However, in 

2002-2004, the blue wing of both lines became steeper than the 
red one and it did not contain any noticeable features. 

2) In 2005 (May-June), when the nucleus of NGC 4151 
was in the minimum activity, the line profiles had a double-peak 
structure with two distinct peaks (bumps) at radial velocities of 
(-2586; +2027) km s-i in H/? and (-1306; +2339) km s"' in Ha 
(see Fig. H 2005_05 and 2005_06). In principle, the two-peak 
structure in the H/3 profiles is also seen in spectra of 1999-2001: 

at y, 1500 km s"' the blue and Vr ~ 500 km s"' the red 

peak. But in this case the blue peak may be caused by a broad 
absorption line at the radial velocity ( 400 km s"'). 

3) In 2006 the line profiles changed dramatically - the blue 
wing became flatter than in previous period, while the red wing 
was very steep without any feature at Vr > 2300 km s"' . 

4) In 2002 a distinct peak (bump) appeared in the red wing 
of the Ha and Hy6 lines at the radial velocity Vr ~ 3000 km s ' . 
The radial velocity of the red peak decreases: in 2002-2003 it 
corresponds to ~ 3100 km s"' and in 2006 to ~ 2100 km s"'. 
This effect is well seen in Fig. ID especially in the Ha line profile. 
Table[T|gives the obtained radial velocities of red peak measured 
in the Ha and H/3 profiles for which the peak was clearly seen. 
Radial velocities of the Ha and H/? red peak, measured in the 
same periods, are similar (the differences are within the error- 
bars of measurements). The mean radial velocity of the red peak 
decreased by ~ 1000 km s"' from 2002 to 2006 (Table [B. It 
is not clear if the red peak is shifting along the line profile or 
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Fig. 6. The absorption seen in years 1996, 1997, 1998 and 2000 
(from top to bottom). The H/3 line after subtraction of the narrow 
lines is shown with solid lines. The vertical line marks the radial 
velocity of -400 km/s. 



it disappears and again appears as a new red peak at another 
velocity. 

3.2. The absorption and emission features in Ha and H/3 line 
profiles 

As mentioned above (§3.1), the Ha and H/3 line profiles have the 
absorption and emission features. Distinguishing between these 
features is very difficult. The question is are the "bumps" an in- 
trinsic property of the broad Ha and Hy6 lines or is it just the 
strengthening of the absorption features? Furthermore, an open 
question is the origin of these features, i.e. is there an intrinsic 
mechanism in the BLR which creates these "bumps"? We should 
mentioned here that in the central part of both lines, the residu- 
als of the narrow components can affect these relatively weak 
absorption/emission features. 

To study it, the residuals between some H/3 broad lines with 
prominent (noticeable) bumps from particular two successive 
months were obtained. An example of these residuals are pre- 
sented in Fig. |5] (available electronically only). The noticeable 
residual bumps (without an absorption-like feature) at Vr from (- 
2000) km s"' to (-1000) km s"' in 1999-2001, and at V,- from 
3500 km s-i to 2500 km s"' in 2002-2006 are seen well. It 
seems that the absorption changes slowly and during several next 



Table 1. The peak shifts in the red wing of the Ho- and H/3 Unes. Columns: 1 - year and month for Ha (e.g., 2002_1_3 - 2002, January 
and March); 2 - y,(Ha), Ha peak velocity in the red wing (km s~'); 3 - year and month for H/?; 4 - Vr(li(3), Hp peak velocity in the 
red wing (km s '). The line after each year (or interval of years) gives the mean yr(Ha) and Vr(li/3) and their standard deviations. 
The last line gives the maximal shift of the red peak (AVj) from 2002 to 2006. 
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months it remains constant. Therefore this absorption disappears 
in profile residuals (see Fig.|5]l. Consequently, it seems that emis- 
sion bumps observed in the H/3 and Ha line profiles (seen in the 
residuals in Fig.|5]l are mostly an intrinsic property of the broad 
emission-line profile. 

The strong absorption features are also present in 1996-2001 
Hj3 spectra. In this period we observed the dip and broad absorb- 
tion at radial velocity that changes from ~-1000 km/s in (1996- 
1998) to —400 km/s in 1999-2000 (see Fig. H available elec- 
tronically only). This velocity corresponds to the minimum of 
the absorption band, but its blue edge extends to higher veloc- 
ities (—1800 km/s in 1998 and —1170 in 1999-2000). Fig. |6] 
shows some observed individual spectra and their broad compo- 
nent where the blue absorbtion is well resolved. The blue shifted 
absorption is probably coming from an outflowing material. 

It is interesting to note that the higher radial velocity ( — 
1000 km/s, observed in 1996-1998) is appearing for a higher 
continuum flux level, while the smaller velocity (~ -400 km/s) 
has been detected when the continuum flux decreases 3-6 times, 
i.e. we confirm the results reported by Hutchings et al. (2002), 
who found the same trend that the outflow velocity increases 
with the continuum flux. 



3.3. Asymmetry of the broad Ha and H/3 line profiles 

We measured the full width at half maximum (FWHM) of the 
broad lines from their month-averaged profiles and determined 
the asymmetry (A) as a ratio of the red and blue parts of FWHM, 
i.e. A=WredAVbiue, where Wred and Wbiue are the red and blue 
half- widths at maximal intensity (see Fig.|7] available electroni- 
cally only) with respect to the position of the narrow component 
of Ha and Hy6. As we mentioned above, there are residuals in 
the center of the Ha and H/3 lines due to residuals from the sub- 
traction of narrow components, that can affect the FWHM and 
A measurements. Therefore we first smoothed the line profiles. 



in order to avoid artificial peaks from the residuals (see Fig.]?}, 
before measuring the FWHM and asymmetry. The two indepen- 
dent measurements of FWHM and A were performed. We deter- 
mined the averaged continuum and its dispersion for each month 
and the averaged Julian date from the spectra which were used to 
construct the month-averaged profiles. The measurements of the 
FWHM, asymmetry and continuum during the whole period of 
monitoring (1996-2006) are presented in Table|2](available elec- 
tronically only) and in Fig. [8] In Table|2]we give averaged values 
of the FWHM and A from two independent measurements and 
their dispersions. 

The FWHM of H/3 was almost always larger than the Ha one 
(see upper panel of Fig.|8]l. The asymmetry of both lines (middle 
panel of Fig.[8]l has been gradually increasing from 1996 to 2006 
and it slightly anticorrelates with the variations of the continuum 
(bottom panel in Fig.|8]l. The largest values of the FWHM and an 
outstanding red asymmetry of both lines (A>1.2) was observed 
in 2002-2006. We calculated average values of the FWHM and 
A for each year, as well as for the whole monitoring period (they 
are given in Table|2li. FWHMs and asymmetries obtained in this 
work from measurements of the month-average profiles are sim- 
ilar to the results given in Paper I and difference between them 
is within the error bars. As one can see from Table|2]the FWHM 
of both lines was varying rather considerably from year to year 
(AFWHM~ 500- 1500km s"' ). These lines were the narrowest 
in 2000-2001 (FWHM ~ 4000 - 5000 km s"') and the broad- 
est in 2005 (~ 7000 kms"') (see Table |2ll. At the same time 
the H/3 FWHM was always broader than Ha by ~ 1000 km s"' 
on average. The asymmetry was varying in different ways: in 
1996-1997 the blue asymmetry was observed in H/3 (A~0.85) 
when the continuum flux was maximum; in 1998-2000 Hp was 
almost symmetric, and from 2001 to 2006 the red asymmetry 
appeared (A>1.2). In 1996-2001, the blue symmetry (A~0.8) 
was observed in Ha, and in 2002-2006 it was a red asymmetry 
(A>1.2). 
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Fig. 10. Light curves of the different Ha (upper panel) and Hfi (bottom panel) line segments. The segments in the blue wing (numbers 
from -5 to -1 from Table O are marked with crosses (x) and in the red wing (numbers from +5 to +1 from Table |3]l with plus (+). 
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comparable with the variation of the central part of the Une, is presented with the solid line. 



Also, we tried to find correlations between the FWHM and 
A with the continuum flux and found that the FWHM practically 
does not correlate with the continuum (r~ 0.0; -0.1). In the case 
of the asymmetry A, there is an indication of anticorrelation, 
but it should be also taken with caution since there is a large 
scatter of points on A vs. continuum flux plane, especially in the 



case of low continuum fluxes < 5.3 x lO^'^^erg cm^^s^'A"' 
when the measured asymmetry reaches the highest values. Note 
here that photoinoization model predicts that the Balmer lines 
should be broader in lower continuum states and narrower in 
the higher continuum states (see Korista & Goad 2004), as due 
to larger response in the line cores one can expect that Balmer 
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Fig. 11. The segment to segment response, where the first period (Period I, 1996-1999) is denoted with open circles, the second 
period (Period II, 2000-2001) with full ones, and the third period (Period III, 2002-2006) with open triangles (Paper I). The flux is 
given in erg cm"^ s"' . 



Table 3. The beginning and ending radial velocities, Vbeg and Vend, in km/s for different segments in the line profiles. 
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lines became narrower in higher continuum states. As one can 
see from Table |2] there is no trend that FWHM is significantly 
narrower in the high continuum state. 

3.4. Light curves of different line segments 

In Paper I we obtained light curves for the integrated flux in lines 
and continuum. To study the BLR in more details, we are starting 
in this paper from the fact that a portion of the broad line profile 
can respond to variations of the continuum in different ways. 
Therefore, we divided the line profiles into 1 1 identical profile 
segments, each of width ~1000 km s"' (see TableO. 

The observational uncertainties were determined for each 
segment of the Hfi and Ha light curves. In the uncertainties we 
include the uncertainties due to: the position angle correction, 
seeing correction procedure and aperture correction. The meth- 
ods to evaluate these uncertainties (error-bars) are given in Paper 
I. The effect of the subtraction of the template spectrum (or the 
narrow components) has been studied by comparing the flux of 
pairs of spectra obtained in the time interval from to 2 days. 
In Table 13] (available electronically only) we presented the year- 
averaged uncertainties (in percent) for each segment of Ha and 
H/3 and mean values for all segments and corresponding mean- 
year flux. To determine the error-bars we used 44 pairs of Ha 
and 68 pairs of H/3. As one can see from Table HI for the far 
wings (segments +5) the error-bars are greater (~10%) in H/3 
than in Ha (~6%). But when comparing the error-bars in the far 



red and blue wings, we find that the error-bars are similar. Also, 
higher error-bars can be seen in the central part of the Ha due 
to the narrow line subtraction. Fig. |9] (available electronically 
only) shows the distributions of the error-bars as a function of 
the line flux for segments ha+5 and haO and ha+ 1 . It can be seen 
that in the case of haO and ha-i-1 there is a slight anticorrelation 
with flux and two points, corresponding to the very small flux 
(F< 4-10"'"'erg/cm"^sec"' in 2005), have the highest error-bar 
of (40-70)%. 

We constructed light curves for each segment of the Ha and 
H/3 lines. Fig.[TO]presents light curves of profile segments in ap- 
proximately identical velocity intervals in the blue and red line 
wings (segments from 1 to 5, where larger number corresponds 
to higher velocity, see Table O and for the central part (0 in 
Table [3] or Ha_c, Hj0_c in Fig. \W[ corresponding to the inter- 
val ~ +500 km s"'). To compare the segment variation with the 
continuum we plot (as a solid line) the continuum flux variation 
into central part (see Fig.fTOb. 

In 1996-1997 the blue segments 2 and 3 were slightly 
brighter than the red ones, while the segments 1, 4 and 5 were 
similar to the red ones. In 1998-2001 the blue segments 4 and 
5 (3500-5500 kms"', regions the closest to the BH) and the 
segment 3 (in 2002-2006) of both lines were essentially fainter 
than the red ones. In 1998-2004 the blue segments 1 and 2 (500 
km/s - 3500 km/s) were close to the corresponding red ones or 
slightly fainter (see Fig.fTOb. 
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Fig. 8. Variations of the FWHM (upper panel), asymmetry (mid- 
dle panel) in Ha (denoted with crosses) and Hj3 (denoted with 
plus) broad lines, and of the continuum flux at /15100A (bot- 
tom panel) in 1996-2006. The abscissae shows the Julian date 
(bottom) and the corresponding year (up). The continuum flux is 
given in lO^'^^erg cm"^ s"'A"'. 



3.5. The line-segment to line-segment flux and 
continuum-line-segment relations 

The line-segment to line-segment flux and continuum-line- 
segment relations for the Ha and H/? lines are practically the 
same, therefore here we present only results for the Ha line. 
First, we are looking for relations between the Ha segments 
which are symmetric with respect to the center (i.e. segments -1, 
and 1; -2 and 2, ... -5 and 5). In Fig.[TT]we present the response 
of symmetrical segments of Ha to each other for the three peri- 
ods given in Paper I. As can be seen from Fig. [TT]the symmetric 
segments are pretty well correlated, with some notable excep- 
tions: a) weaker response of the red to the blue wing in the II and 
III period and b) the apparent bifurcation in the Ha2/Ha-2 and 
Ha3/Ha-3 plots. This appears to be associated with Period III, re- 
lated to the appearance of the 3000 km/s "red bump". This also 
supports that lines are probably formed in a multi-component 
BLR and that the geometry of the BLR is changing during the 
monitoring period. 

Additionally, in Fig. [12] we present the response of dif- 
ferent Ha segments to the continuum flux. As it can be seen 
the responses are different: in the far wings, the response to 
the continuum is almost linear for the red wing (segments 4 
and 5) and for a fraction of the blue wing (segment -4), but 
for the far blue wing (-5500 to -4500 kms"') there is prac- 
tically no response for < 7 x lO^'^^erg cm^^s^'A"'. For 



higher continuum fluxes, the far blue wing has a higher flux, 
but it seems also that there is no linear relation between the 
line wing and the continuum flux. On the other hand, the cen- 
tral segments (from - 3500 to 3500 km s"') have a similar re- 
sponse to the continuum like the Hfi and Ha total line fluxes 
(see Paper I) - a linear response for the low continuum flux 

< 7 X 10 '''erg cm^^s ' A ') and no linear response for the 
high continuum flux > 1 x lO^'^'erg cm^^s^'A"') The linear 
response indicates that this part of line (red wings and central 
segments at the low continuum flux) is coming from the part of 
the BLR ionized by the AGN source, while the blue and partly 
central part of the line could partly originating from a substruc- 
ture out of this BLR (and probably not photoionized). Note here 
that the photoionization in the case of a mix of a thin (fully ion- 
ized) and thick (ionization-bounded) clouds can explain the ob- 
served non-linear response of emission lines to the continuum in 
variable as seen in the central parts of broad lines (see Shields et 
al. 1995), i.e. in the case of the optically thick BLR the detailed 
photoionization models show that the response of the Balmer 
lines declines as the continuum flux increases (see Goad at al. 
2004, Korista & Goad 2004). In our case we found that the flux 
in wings (except the far blue wing) have almost linear response 
to the continuum flux, while the central parts show non-linear 
response to the higher continuum flux. Also the response to the 
continuum flux of the far blue wing (-5) and far red wing (-1-5) 
is very different. It may indicate different physical conditions in 
sub-regions or across the BLR. 

In Fig. [13] (available electronically only) the flux of differ- 
ent Ha segments as a function of the flux of the central segment 
(HaO) are presented. It can be seen that there are different rela- 
tions between different Ha segments and HaO: for segments near 
the center (Ha-1, Ha-2, Hal and Ha2) the relation is almost 
linear, indicating that the core of the line originates in the same 
substructures (Fig.fTTt: segments in the Ha wings (the near blue 
wing Ha-3 and the red wing Ha3, Ha4 and Ha5) also show a lin- 
ear response to the central segment (but the scatter of the points 
is larger than in previous case) which also indicates that a por- 
tion of the emission in the center and in these segments is com- 
ing from the same emission region. On the contrary, the far blue 
wing (Ha-4 and Ha-5) responds weakly to the line center, espe- 
cially Ha-5 that shows practically independent variations with 
respect to the central segment. 

4. The Balmer decrement 

4.1. Integral Balmer decrement 

From month-averaged profiles of the broad Ha and li/3 lines 
we determined their integrated fluxes in the range between - 
and +5500 km s ' of radial velocities. We call "Integral Balmer 
decrement" (BD) the integrated flux ratio F(Ha)/F(li/3). 

Fig. [14] shows the behavior of the integrated BD (upper 
panel) and continuum flux at the wavelength 5 100 A (bottom 
panel). In 1999-2006 an anticorrelation between the changes of 
the integrated BD and continuum flux was observed. It was es- 
pecially noticeable in 1999-2001. 

Table |5]presents year-averaged values of the BD and contin- 
uum flux determined from month-averaged profiles in each year 

We found that in 1996-1998 the continuum flux was rather 
large and was changing within the limits F^ ~ (6 - 12) x 
10"'"' erg cm^^s^'A"', but the BD was practically not chang- 
ing. In Paper I we akeady noted the absence of correlation be- 
tween the continuum flux and the integrated flux of the broad 
lines for the above-mentioned flux values. Also, we found (see 
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Fig. 12. The flux of diff'erent segments of the Hne as a function of the continuum flux. The Hne flux is given in erg cm s and the 
continuum flux in erg cm^^s^'A"'. 



Figure fl4l that the Balmer decrement is systematically higher in 
1999-2001. 



Table|5]also gives average data of the BD and continuum flux 
for the periods of 1996-1998, 1999, 2000-2001, 2002-2006. It 
is evident that in the period of 1996-1998 the BD did not vary 
within the error bars in spite of strong variations of the contin- 
uum. Radical changes (the increasing of the BD) started in 1999. 



The BD reached its maximum in 2000-2001. Then in 2002- 
2006 average values of the BD coincide with those of 1996- 
1998. Thus, it can be concluded with confidence that from 1999 
to 2001 we observed an obvious increase of the BD. 



4.2. Balmer decrement of different profile segments 

From month-averaged profiles of Ha and we determined the 
BD for the segments of Table|3] Fig.[T5](available electronically 
only) shows the BD variation in a 2D plane (year-V^), and Fig. 
[T6]gives the changes of the BD of different segments of the pro- 
files during the whole monitoring period. 

It can be seen from Figs. [B] and \T6\ that in 2000-2001 
(JD=5 1200-52300) the values of the BD in diff'erent segments 
were always on average noticeably larger than in the other years, 
and the blue wing had always on average a noticeably larger BD 
than the red one. In 1996-1998 BDs of the blue and red wings 
practically coincided, and in 2002-2006 they either coincided or 
the BD of the blue wing was slightly larger. 
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Fig. 14. Variations of the integrated B aimer decrement 
BD-F{Ha)fF(H/3 (upper panel) and of the continuum flux at 
/15100A (bottom panel) in 1996-2006. The abscissae gives the 
Julian date (bottom) and the corresponding year (top). The con- 
tinuum flux is in units 10^'^ergcm^^s"' A"'. The vertical lines 
coiTespond to yeai^s 1999 and 2002. 
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Fig. 16. Variations of the BD of dififerent segments of the line 
profiles and of the continuum flux (bottom panel) in 1996-2006. 
The BD of segments in the blue wing (numbers from -5 to -1 
from Table |3]l are denoted with plus (+), in red wing (numbers 
from -1-5 to -i-l from Table 13) with crosses (x), and of central 
part (number 0) with open circles. The vertical lines correspond 
to years 1999 and 2002. The abscissae shows the Julian date 
(bottom) and the corresponding year (top). The ordinate shows 
the BD for different segments of the line profile. The continuum 
flux is in units 10"'^ erg cm"^ s"' A"' . 



During the whole monitoring period the value of the BD 
along the line profile in central segments (in segments 0,+l cor- 
responding to the velocity range from -1500 km s ' to 1500 
km s~') were considerably larger (by ~ 1.5) on average than 
at the periphery (segments from +2 to +5 corresponding to 



the velocity range from ((-5500)-(-1500) kms"') to (-H5500)- 
(-Hl500)kms-i)). 

4.2.1 . Balmer decrement variations as a function of tine radial 
velocity 

As an illustration of the variation of the Balmer decrement with 
the radial velocity during the monitoring period, we show in Fig. 
[17] the BD as a function of the radial velocity. Since the BD 
vs. velocity remains the same during one year, we give a few 
examples just for month-averaged spectra for each year of the 
monitoring period. 

On the other hand, the behavior of the Balmer decrement 
as a function of the radial velocity differs in different years. 
As a rule, in 1996-1998 the maximum of the BD was at ~ 
-1000 km s"', while the BD was slowly decreasing in the ve- 
locity range from to 1000 km s"', and sharply decreasing in 
the region > +1000 kms"', usually more strongly in the re- 
gion of negative velocity. Recall here the results obtained from 
photoionization model by Korista & Goad (2004) where they 
found velocity dependent variations in the Balmer decrement. 
They found that the Balmer decrement is steeper in the line core 
than in line wings, as we obtained in some cases (see Fig. [TtI i. 
but it is interesting that in all periods this peak is offset from the 
central part to the blue side and also there are several cases where 
the Balmer decrement is steeper in the wings. On the other hand 
in the case if velocity field is dominated by the central massive 
object, one can expect symmetrical BD in the blue and red part 
of velocity field (Korista & Goad 2004). In our case we obtained 
different asymmetrical shapes of BD vs. velocity field. The BD 
seems to have a systematic change in behavior starting around 
2002 - i.e. corresponding to period III, when the "red bump" ap- 
pears, showing a two maxima in the BD vs. velocity field. This 
may indicate that velocity field in the BLR is not dominated by 
central massive black hole, i.e. it is in favor of some kind of 
streams of the emitting material as e.g. outflow or inflow. 

In 1999-2001 the maximum values of the BD ~ (6-8) were 
observed in the velocity region +1000 km s"', while a steeper 
decrease of the BD was more often observed in the region of the 
positive velocity. The change of the BD relative the radial ve- 
locity in 2002-2006 differed strongly from those in 1996-2001. 
In these years 2 peaks (bumps) were observed in the BD distri- 
bution: 1) at radial velocities from -2000 to -1000 with larger 
values of the BD; 2) at +3000 km s"' with somewhat smaller 
(by ~ (0.5 - 1.0)) values of the BD. 

4.3. Balmer decrement and helium line ratio 

In order to probe the physical conditions in the BLR, we studied 
the flux ratio of He II /i4686 and Hel /i5876 broad lines. From 
all the available spectra, we selected only 21 spectra where the 
broad helium lines could be precisely measured, and where the 
two helium lines were observed on the same night. Note that in 
the case of the minimum activity, the broad component of the He 
II /14686 line could not be detected at all. 

We use the helium lines He II /14686 and He I /15876 since 
these two lines come from two different ionization stages and, 
thus, are very sensitive to changes of the electron temperature 
and density of the emitting region (Griem 1997). 

In Fig. [18] we give the Hell/Hel vs. continuum flux (first 
panel), BD vs. Hell/Hel (second panel) and He II/Hy6 vs. contin- 
uum flux (third panel). As it can be seen from Fig.[T8](first panel) 
there is a good correlation between the continuum flux and the 
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Fig. 17. Variations of the Balmer decrement BD=F(Ha)/F(I-lj6) as a function of the radial velocity for month-averaged spectra each 
year of the monitoring period. The abscissae gives the radial velocities relative to the narrow components. The month-averaged 
Juhan date and the corresponding year are given at the top of each plot. 



Fig. 18. Left: The variation of the ratio of the heHum He II /14686 and He I /I5876 lines as a function of the continuum flux. Middle: 
The BD as a function of the helium line ratio. Right: The variation of the ratio of the helium He II /14686 and H/3 lines as a function 
of the continuum flux. The continuum flux is in units 10 '"^erg cm^^sec ' A 



helium line ratio (the correlation coefficient r=0.81) which in- 
dicates that these lines are probably coming from the region 
photoionized by the continuum source. On the other hand, the 
Balmer decrement decreases as the ratio of the two helium lines 
increases (the anti-correlation is a bit weaker here, the correla- 
tion coefficient is r=0.75), indicating that as the ionization gets 
stronger, the BD decreases (Fig. [18] second panel). 

The photoionization model (Korista & Goad 2004) shows 
that the Balmer lines should show less flux variation with con- 
tinuum state than He I, that is less varying than He II. In this 
case one can expect that He II/Hy6 changes much more than the 
He II/He I. The observed changes in the He II/H/? flux ratio (See 
Fig. [18] third panel) is probably not due to differences in the 
way these two lines respond to the level of the ionizing contin- 
uum flux (see Shapovalova et al. 2008, Fig. 10 in the paper), but 
rather to changes in the shape of the ionizing SED. In our case, 
this ratio drops about 4 times from the high to low continuum 
states, while the He II/He I line ratio changes around 3 times. 

There is some connection between the helium and Balmer 
lines (Fig. [18] left), but the physical properties of the emitting 
regions of these lines are probably different. One of the explana- 
tion of this correlation may be that a part of the Balmer lines is 
coming from the same region as He I and He II lines. Since the 
ratio of He II/ He I is sensitive to the change in temperature and 
electron density, in the case of higher He II/ He I ratio the ion- 
ization is higher and population of higher level in Hydrogen has 
higher probability, consequently the ratio of HafH/3 is smaller 
than in the case of lower He II/ He I ratio. 



5. Summary of the results 

In this paper we investigated different aspects of broad Ha and 
H/3 line variations (shapes, ratios, widths and asymmetries), so 
in this section we outline the summaries of our obtained results. 

5.1. Summary of the line profile variations 

i) The line profile of Ha and H/3 were changing during the 
monitoring period, showing blue (1996-1999) and red (2002- 
2006) asymmetry. We observe bumps (shelves or an additional 
emission) in the blue region in 1996-1997 at Vr ~ -4000 and 
-2000 km s-\ in 1998-99 at V, ~ -2000 km s"' , and in 2000- 
2001 at -2000 and -500 km s"' akeady. However, in 2002-2004 
these details disappeared in the blue wing of both lines. It be- 
came steeper than the red one and did not contain any noticeable 
features. In 2002 a distinct bump appeared in the red wing of 
both lines at -1-3100 km s"'. The radial velocity of a bump in the 
red wing changed from H-3100 km s ' in 2002 to ~ 2100 km s ' 



Table 5. Year-averaged and period-averaged variations of the in- 
tegral Balmer decrement (BD) and continuum fluxes. Columns: 
1-year or year-intervals; 2- F(Ha)/F(Hy6), Ha and Hfi flux ra- 
tio or integrated Balmer Decrement (BD); 3 - cr(F(Ha)/F(Hy6)), 
the estimated F(Ha)/F(Hy6) error; 4 - F(5100), continuum flux 
at /15100 A; 5 - cr, estimated continuum flux error. The last four 
lines give the mean BD, mean continuum flux and their esti- 
mated errors in year-periods. 



year 


¥(lla)/F{Hj3)±(T 


F(5100)+cr;i 


1996 


3.65 


0.51 


11.18 


1.55 


1997 


3.39 


0.36 


7.62 


1.03 


1998 


3.58 


0.35 


6.08 


1.19 


1999 


4.49 


0.74 


5.12 


1.62 


2000 


5.01 


0.62 


2.75 


1.33 


2001 


5.43 


1.07 


1.87 


0.47 


2002 


3.67 


0.44 


3.59 


0.29 


2003 


3.16 


0.27 


4.84 


0.58 


2004 


3.70 


0.29 


2.44 


0.57 


2005 


3.96 


0.54 


1.99 


0.18 


2006 


3.58 


0.36 


2.60 


0.80 


mean (1996-1998): 


3.54 


0.13 


8.29 


2.62 


mean 1999(Jan.-Apr.): 


4.21 


0.59 


5.89 


0.56 


1999 (Dec): 


5.33 




2.78 


0.12 


mean (2000-2001): 


5.22 


0.30 


2.31 


0.62 


mean (2002-2006): 


3.61 


0.29 


3.09 


1.14 



" in units of 10 '''erg cm -s 'A ' 



in 2006. Possibly, the appearance of the blue and red bumps are 
related with a jet-component. 

ii) In 2005 (May- June), when the nucleus of NGC 4151 was 
in the minimum activity state, the line profiles had a double- 
peaked structure with two distinct peaks (bumps) at radial veloc- 
ities of -2586; +2027 km s"' in Hfi and -1306; -h2339 km s"' 
in Ha. 

iii) In 1996-2001 we observed a broad deep absorption line 
in Hy6 at a radial velocity that took values from V,- ~ -1000 to 
y, ~ -400kms-i. 

iv) the FWHM is not correlated with the continuum flux, 
while the asymmetry tends to anti-correlate with it (coefficient 
correlation r — 0.5). The FWHM and asymmetry in li/3 are larger 
than in Ha. It may be explained with the fact that H/3 originates 
closer to the SMBH, and thus has larger widths. 

v) We divided the line profiles into 1 1 identical segments, 
each of width 1000 km s"' and investigated the correlations be- 
tween segments vs. continuum flux and segments vs. segments 
flux. We found that the far red wings (from 3500 to 5500 km s"') 
and central (at Vr - ±3500 km s"') segments for the continuum 
flux Fc <1 X lO^'^^ergcm^^s"' A"' respond to the continuum al- 



most linearly, and this segments are probably coming from the 
part of the BLR ionized by the AGN source. 

vi) The central (+3500 km s"') segments for > 7 x 
10"''*ergcm"^s 'a ' do not show any linear relation with the 
continuum flux. Probably in periods of high activity this seg- 
ments of line is partly originating in substructures which are not 
photoioninized by the AGN continuum. 

vii) The far blue wing (from -5500 to -4500) seems to origi- 
nate in a separate region, so it does not respond to the continuum 
flux variation as well as to the variation of other line segments, 
except in the case of high line flux where it responds to the far 
red wing (see Fig.fTZt. This may indicate that the far blue wing is 
emitted by a distinct region which is not photoionized, and also 
that the emission is highly blueshifted (as e.g. an outflow with 
velocities >3500 kms"'). 

viii) The far red wing is very sensitive to the continuum flux 
variation, and is thus probably coming from the part the closest 
to the center, i.e. this part of the line seems to be purely pho- 
toionized by the AGN source. 

From all the facts mentioned above, one can conclude that 
the broad lines of NGC 4151 are produced at least in three kine- 
matically (and physically) distinct regions. 

5.2. Summary of the BD variations 

The Balmer decrement (BD=F(Hff)/F(HjS)) varied from 2 to 8 
during the monitoring period. It is interesting that BDs are quite 
different along line profiles as well. In 1996-1998 there was no 
significant correlation between the BD and the continuum flux 
(the continuum flux was ~ (6 - 12) x lO^'^^erg cm^^s"' A"'). 
In 1999-2001 maximal variations of the BD were observed, es- 
pecially in the blue part of lines. The maximal value of the 
BD along the line profiles strongly differed: in 1996-2001 the 
BD was maximum at Vr ~ +1000 km s"^ and in 2002-2006 
the BD had 2 peaks - at velocity from -2000 to -1000 and at 
Vr ~ 3000 km s"' with somewhat smaller (by ~(05-l)) values of 
the BD. In the last case it is possible that the second bump (the 
fainter one) is caused by the interaction between the receding 
sub-parsec jet and environment. 

The different values of the BD observed during the monitor- 
ing period (as well as different values of the BD along the pro- 
files) also indicate a multicomponent origin of the broad lines. 
Such different ratios may be caused by absorption, but also by 
different physical conditions in different parts of the BLR. 



6. Discussion 

It is interesting to compare our results with those found in the 
UV and X-radiation. Crenshaw & Kraemer (2007) found a width 
of 1 170 km s"' (FWHM) for the UV emission lines significantly 
smaller than our results for the Balmer lines (FWHM~ 5000- 
7000 km s around 5-6 times smaller). This can be interpreted 
as the existence of an intermediate component between the broad 
and narrow (emission) line regions (see e.g. Popovic et al. 2009). 
In the same paper, they found an evidence that the UV emission 
lines arise from the same gas responsible for most of the UV 
and X-ray absorption. The absorption can be seen in outflow at 
a distance of 0.1 pc from the central nucleus. Also, it is interest- 
ing that Crenshaw & Kraemer (2007) favor a magnetocentrifugal 
acceleration (e.g., in an accretion disk wind) over those that rely 
on radiation or thermal expansion as the principal driving mech- 
anism for the mass outflow. Obscuration should play an impor- 
tant role, but we found that while an absorption can clearly be 



detected in the HjS line, the Ha line displays a small amount of 
absorption (see Figs. |2]and[3]l. Absorption in the H/3 line was 
previously reported by Hutchings et al. (2002), therefore obscu- 
ration of the optical continuum should be present in NGC 4151. 
The location of the obscuring material was estimated by Puccetti 
et al. (2007). They analyzed the X-ray variability of the nucleus. 
They found that the location of the obscuring matter is within 
3.4x10^* Schwarzschild radii from the central X-ray source and 
suggested that absorption variability plays a crucial role in the 
observed flux variability of NGC 4151. 

Let us discuss some possible scenarios for the BLR. As we 
mentioned above, the outflow is probably induced by the magne- 
tocentrifugal acceleration, and starts very close to the black hole. 
If we adopt a black hole mass around 4x10^ Mq (obtained from 
stellar dynamical measurements, see Onken et al. 2008, which 
is in agreement with Bentz et al. 2006), the acceleration (and 
line emission) is starting at ~10^^ pc from the black hole and 
the outflow is emitting until ~ 0.01 pc; taking into account the 
absorption velocities (e.g. Hutchings et al. 2002) around -1000 
km s"'. 

As it was mentioned above, the broad Ho- and H/5 lines show 
different widths and asymmetries during the monitoring period, 
indicating a complex structure of the BLR. Also, we found that 
the line profiles of Ha and H/3 could be different at the same 
epoch. Therefore one should consider a multi-component origin 
of these lines. 

To propose a BLR model, one should take into account that: 
a) there is an absorption component in the blue part of lines indi- 
cating some kind of outflow that may start in the BLR (see §3.2). 
Also Crenshaw & Kraemer (2007) reported a mass outflow from 
the nucleus of the NGC 4151, but they confirmed the observed 
outflow in the Intermediate Line Region (ILR).; b) the flux in the 
far red wing correlates well with the continuum flux, indicating 
that it originates in an emission region very close to the contin- 
uum source, i.e. to the central black hole; c) the flux of the far 
blue wing does not correlate with the continuum, that may indi- 
cate some kind of shock wave contribution to the Balmer line. 

All these facts indicates that an outflow probably exists in 
the BLR. In this case a complex line profiles (with different fea- 
tures) can be expected due to changes in the outflowing structure, 
as it is often seen in the narrow lines observed in the jet induced 
shocks (see e.g. Whittle & Wilson 2004). Of course we cannot 
exclude that there is the contribution of some different regions to 
the composite line profile, as e.g. there can be also contribution 
of the ILR which may be with an outflow geometry. In forth- 
coming paper we will try to adjust the observational results with 
those predicted by various models of the kinematics and struc- 
ture of the BLR (a bi-conical outflow, an accelerating outflow 
(wind), a Keplerian disk, jets, etc). 

It is interesting to note that the line profiles are changing 
during monitoring period, especially after 2002, when the 'red 
bump' appears. After that the asymmetry of both lines (see Table 
2) and BD show different behavior than in the period 1996 - 
2002. Moreover, in the III period, the line profiles of Ho- and 
H/3 very much changed in the shape from double peaked pro- 
files (see Figs. 1 and 2) to a quite asymmetrical profiles (as was 
observed in 2006). 

The integrated Balmer decrement was maximum in 1999- 
2001 (see Fig. fl4l i. The BD changes shape from 2002, show- 
ing two peaks in the BD vs. velocity field profile. Most prob- 
ably, this is connected with strong inhomogeneities in distri- 
bution of the absorbing material during different periods of 
monitoring. As the line flux in Ha and H/3 at small fluxes 
(< 7 X lO^'^^erg cm^^s^'A"') correlates well with that of the 



continuum (Paper 1), we infer that the change of the integrated 
Balmer decrement in 1999-2001 is also caused, at least partially, 
by changes in the continuum flux. Indeed, when the ionization 
parameter decreases for a constant density plasma, an increase of 
F(HQr)/F(HyS) intensity ratio is expected (see for instance Wills et 
al. 1985). This is due to the decrease of the excitation state of the 
ionized gas: the temperature of the ionized zone being smaller, 
the population of the upper levels with respect to the lower ones 
decreases. In 1996-1998 BDs did not correlate with continuum, 
i.e. in this case the main cause of BD variations is not related 
to the active nucleus and probably the shock initiated emission 
is dominant. We detected that the FWHM of lines does not cor- 
relate with continuum. This confirms our assumption that broad 
Unes can formed in several (three) different subsystems or that 
the emission is affected by an outflow that produces shock initi- 
ated emission. 



7. Conclusions 

This work is subsequent to the Paper I (Shapovalova et al. 2008) 
and is dedicated to a detailed analysis of the broad Ha and H/3 
line profile variations during the 11 -year period of monitoring 
(1996-2006). From this study (Section 3) it follows that the BLR 
in NGC 4151 is complex, and that broad emission hues are the 
sums of components formed in different subsystems: 

1) the first component is photoionized by the AGN contin- 
uum (far red line wings, Vr ~ +4000 and +5000 km s"' and 
central (at Vr - +3500 km s ') segments for continuum flux 
Fc < 7 X lO'^'^ergcm'^s"' A"'). This region is the closest to the 
SMBH. 

2) the second component is independent of changes of the 
AGN continuum (far blue line wings, Vr ~ -5000 km s"'; 
-4000 km s"^). It is possibly generated by shocks initiated by 
an outflow. 

3) the third component, where the central parts of lines 
(Vr < 4000 km s"') are formed, in high fluxes Fc > 7 x 
10~'"*erg cm"^s"'A ' is also independent of the AGN continuum 
(possibly, outflow and jet). 

Finally, we can conclude that the BLR of NGC 4151 may not 
be purely photoionized, i.e. besides photoionization, there could 
be some internal shocks contributing to the broad lines see also 
discussion in Paper I). At least there are three different subre- 
gions with different velocity field and probably diflFerent phys- 
ical conditions, that produce complex variability in the broad 
lines of NGC 4151 and changes in the line profile (very often 
temporarily bumps). Our investigations indicate that the rever- 
beration might not be vaUd as a tool to determine the BLR size 
in this AGN and that this AGN is not perfect for this method. 
Consequently, the results for the BLR size of NGC 4151 (given 
in the introduction) should be taken with caution. 
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Fig. 1. The decomposition of (upper panel) and Ha (bottom panel) with Gaussians (below) in order to construct narrow line 
templates. The dashed Gaussians correspond to the narrow components taken for the narrow Une templates, while the soUd ones 
correspond to the broad components (down). The observed spectra and rms (up) are denoted with soUd hne. 
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Fig. 2. The month-averaged profiles of the Ha and Hfi broad emission lines in the period 1996-2006. The abscissae shows the radial 
velocities relative to the narrow component of Ha or Hj8. The ordinate shows the flux in units of 10"^'' erg cm"^ s~^ . 
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;. 2. Continued. 
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Fig. 5. The Hy6 broad emission lines of NGC 4151 in particular two sequent months where bumps were strong (red and blue line). 
Their residual is also given (black Une below). 
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Fig. 7. An example of the FWHM and asymmetry measurements. The observed spectra is denoted with the dashed line and the 
smoothed spectra with solid line. 



Table 2. The FWHM and asymmetry of the Ha and H/? broad emission lines in the period of 1996-2006. Columns: 1 - year and month (e.g., 
1996.01 = 1996 January); 2 - Julian date (JD); 3 - the FWHM of HyS; 4 - the asymmetry of 1^8; 5 - the FWHM of Ha ; 6 - the asymmetry of Hor; 
7 - F(5100), the continuum flux at /15100 A in units of 10'''' erg cm"^ s"' A"' and cr(cnt), the estimated continuum flux error in the same units. 



yearjnonth 


JD 


FWHM Bfi + o- 


AHS±cr 


FWHM Ha ± cr 


AHa±cr 


F(5100)±cr 




2400000-1- 


km s"' 




km s~' 




10-^'^ erg s-^A"' 


1 


2 


3 


4 


5 


6 


7 


1996b_01 


50096.0 


6860 ± 566 


0.904 ± 0.043 


5332 ± 324 


0.799 ± 0.030 


9.031 ± 0.322 


1996b_03 


50163.4 


6182 ± 303 


0.877 ± 0.069 


5036 ± 291 


0.782 ± 0.021 


11.219 ± 0.617 


1996b_04 


50200.8 


6183 ± 479 


0.877 ± 0.072 






9.226 ± 0.527 


1996b.06 


50249.3 


6028 ± 434 


0.797 ± 0.060 








12.641 - 


1996b.07 


50277.8 


6367 ± 479 


0.783 ± 0.048 


5058 ± 259 


0.762 ±0.011 


11.822 ±0.473 


1996b.ll 


50402.6 


5659 ± 348 


0.896 ± 0.035 








10.649- 


1997b.03 


50510.9 


5568 ± 305 


0.884 ± 0.048 


4808 ± 420 


0.853 ± 0.019 


8.350 ±0.716 


1997b.04 


50547.9 


5567 ±218 


0.946 ± 0.046 


4512 ± 323 


0.854 ± 0.005 


6.897 ± 0.375 


1996-1997 




6051 ±449 


0.871 ± 0.054 


4959 ± 307 


0.810 ±0.042 


9.979±1.950 


1998b_01 


50838.0 


5844 ± 173 


1.159 ± 0.005 


4831 ± 323 


0.947 ± 0.050 


6.280 ± 0.277 


1998b.02 


50867.4 


6059 ± 130 


1.189 ±0.021 






6.219- 


1998r_04 


50934.5 






5058 ± 259 


0.981 ±0.051 


6.387 - 


1998b.05 


50940.4 


6489 ± 42 


1.131 ±0.016 


5150 ± 388 


0.933 ±0.071 


6.049 ± 0.385 


1998b.06 


50988.3 


6458 ± 86 


1.059 ±0.029 


4922 ± 323 


0.931 ±0.073 


7.974 ± 0.776 


1998b.07 


51025.3 


6214 ± 88 


1.082 ±0.000 


4762 ± 549 


0.935 ± 0.067 


6.205 - 


1998_11-12 


51148.6 


6275 ± 1 


1.062 ±0.028 


4694 ± 387 


0.909 ± 0.053 


4.975 ± 1.011 


1998 




6223 ± 245 


1.114 ±0.054 


4903 ± 176 


0.939 ± 0.024 


6.298±0.880 


1999b_01 


51202.6 


5967 ± 174 


1.063 ± 0.001 


4831 ± 194 


0.863 ± 0.079 


6.245 ± 0.243 


1999b.02 


51223.1 


5721 ± 174 


0.957 ±0.031 






5.873 ± 0.258 


1999b.03 


51260.5 


5475 ± 261 


1.045 ± 0.064 


4398 ± 291 


0.910 ± 0.045 


5.793 ± 0.296 


1999b.04 


51281.5 


5567 ± 218 


1.011 ±0.078 


4010 ± 388 


0.959 ± 0.001 


5.397 ± 0.146 


1999b.06 


51346.4 


5352 + 522 


1.173 ±0.058 






2.753 - 


1999b.l2 


51515.6 


5106 ± 609 


1.023 ±0.032 


4352 ± 226 


0.836 ± 0.045 


2.782 ±0.016 


1999 




5531 ± 298 


1.045 ±0.072 


4903 ± 337 


0.892 ± 0.054 


4.807± 1.603 


2000b.01 


51553.6 






4010 ± 388 


0.744 ± 0.041 


4.155- 


2000.1.2 


51571.6 


4890 ± 566 


0.960 ±0.124 






3.868 ± 0.406 


2000b.02 


51589.5 


4767 ± 566 


1.036 ± 0.090 


4056 ± 194 


0.801 ± 0.055 


3.581 ± 0.029 


2000b.06 


51711.3 






6244 ± 775 


1.014± 0.162 


1.434 ±0.115 


2000.6.7 


51728.6 


6336 ± 262 


1.423 ±0.060 






1.974 ±0.764 


2000b.ll 


51874.1 


3968 ± 739 


1.014 ±0.020 








1.778 ±0.177 


2000b.l2 


51883.6 


3383 ±783 


1.353 ±0.096 


2985 ± 871 


0.805 ±0.013 


1.427 ±0.09 


2001b.02 


51947.9 


6367 ±914 


1.523 ±0.015 


4603 ± 646 


0.907 ±0.106 


1.538 ±0.320 


2001b.ll 


52238.0 


6152 ±261 


1.127 ±0.027 






3.331 ±0.007 


2000-2001 




5123 ± 1199 


1.205 ± 0.224 


4380 ± 1195 


0.854 ±0.107 


2.565±1.193 


2002b.01 


52299.4 


6183 ± 217 


1.284 ± 0.006 


4785 ± 903 


0.957 ±0.134 


3.680- 


2002r_03 


52345.8 






4603 ± 775 


1.173 ±0.124 


3.348 ± 0.862 


2002.3.4 


52357.3 


6060 ± 479 


1.706 ±0.239 






3.645 ± 0.420 


2002r_04 


52368.7 








5150 ± 838 


1.196 ±0.076 


3.942- 


2002b.05 


52398.7 


6336 ± 173 


1.369 ±0.050 


5605 ± 581 


1.203 ±0.005 


3.283 ± 0.021 


2002b.06 


52439.0 


7012 ±521 


1.285 ±0.088 


5582 ± 226 


1.414 ±0.115 


3.903 ± 0.292 


2002b.l2 


52621.0 


6121 ± 305 


1.489 ±0.051 


5970 ± 387 


1.374 ±0.129 


3.467 ±0.157 


2002 




6490 ± 465 


1.427 ±0.177 


5283 ± 528 


1.220 ±0.163 


3.629±0.243 


2003b.01 


52665.9 


5905 ± 348 


1.710 ± 0.112 


5651 ± 129 


1.373 ±0.117 


5.072- 


2003b.03 


52723.8 


6336 ± 434 


1.344 ±0.066 






4.541 - 


2003b.05 


52777.1 


5906 ± 87 


1.559 ±0.009 


5287 ± 324 


1.260 ±0.027 


5.283 ± 0.205 


2003b.06 


52813.7 


5660± 349 


1.630 ±0.050 


5013 ±388 


1.163 ±0.038 


4.184- 


2003b.ll 


52966.8 


5505 ± 217 


1.358 ±0.082 






3.478 ± 0.193 


2003b.l2 


52995.6 


5414 ± 523 


1.352 ±0.083 


— 


— 


3.515- 


2003 




5788 ± 336 


1.492 ±0.162 


5317 ±320 


1.265 ±0.105 


4.346±0.763 


2004b.01 


53019.5 


5351 ±261 


1.456 ±0.124 


__ 


__ 


3.851- 


2004b.l2 


53349.0 


6889 ± 173 


1.240 ±0.057 


5879 ± 389 


1.422 ±0.076 


2.932 ±0.105 


2004 




6120 ± 1088 


1.348 ±0.153 


5879 ± 389 


1.422 ±0.076 


3.392±0.650 


2005.1.2 


53402.8 


7013 ±783 


1.332 ±0.076 






2.218 ±0.308 


2005b.02 


53417.5 






6586 ± 356 


1.359 ±0.072 


2.040 - 


2005b.05 


53505.3 


8305 ± 262 


1.093 ± 0.043 


6882 ± 582 


1.265 ±0.105 


2.140 ±0.037 


2005b.06 


53535.0 


8427 ± 348 


0.957 ± 0.041 


7635 ± 357 


1.100 ±0.058 


1.841 ± 0.142 



Table 2. Continued. 



yearjnonth 


JD 

2400000+ 


FWHMI^±o- 
km s~' 


AI^±o- 


FWHMHa±£r 
km s"' 


AHar±(r 


F(5100)±o- 
10"" erg s-'A"! 


1 


2 


3 


4 


5 


6 


7 


2005b.l 1 
2005b.l2 


53704.0 
53711.5 


6059 + 130 
5967 ± 88 


1.402 ±0.031 
1.395 ± 0.007 


— 


— 


3.567 - 
2.911 ±0.177 


2005 




7154+1180 


1.236 ±0.200 


7034 ± 540 


1.241 ±0.131 


2.953±0.656 


2006b.01 
2006b.02 
2006b.03 
2006b.04 


53762.4 
53788.0 
53816.9 
53845.1 


6797 ± 218 
6582 ± 174 
6704 ± 86 
5691 ± 305 


1.188 ±0.021 
1.161 ±0.027 
1.117 ±0.002 
1.373 ±0.044 


6015 ±710 
5925 ±518 
5446 ± 613 
5401 ± 356 


1.449 ±0.202 
1.403 ±0.145 
1.546 ±0.222 
1.480 ±0.069 


1.939 ±0.179 
1.926 ±0.016 
2.894 ± 0.293 
3.390 ±0.103 


2006 




6444 ± 509 


1.210 ± 0.113 


5697 ± 318 


1.470 ± 0.060 


2.52>1±0.121 


mean 96-06 




6021 ±851 


1.187 ± 0.237 


5164 ± 873 


1.072 ± 0.245 




Paper I 




6110 ±440 


1.056 ±0.018 


4650 ± 420 


1.00 ±0.023 
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Fig. 9. The flux error measurements (in percent of part line flux) against flux for -5, -1, 0, 1,5 segments of the Ha line. The flux is 
in units of 10 '^ergcm^^sec^'. 



Table 4. The errors of measurements (e+cr) for all line segments (see Table [Sjl of Ha and HyS given in percents. Also, for each 
segment the mean-year flux is given in units lO"'"' erg cm^^s"' . 



Year Ha(-5) Ha(+5) Ha(-4) Ha(+4) Ha(-3) Ha(+3) 





Flux 


(e+cr) 


Flux 


(e+cr) 


Flux 


(e+cr) 


Flux 


(e+cr) 


Flux 


(e+cr) 


Flux 


(e+cr) 


1996 


6.327 


9.6+ 9.7 


7.052 


13.3+10.5 


12.456 


9.9+5.9 


12.494 


11.7+7.1 


22.736 


9.1+6.8 


17.722 


10.9+7.1 


1997 


3.904 


11. 1± 6.1 


4.950 


9.7+ 0.8 


8.100 


6.7+8.8 


10.048 


9.7+3.8 


18.741 


6.7+9.2 


16.634 


9.1+5.3 


1998 


2.497 


11.7± 4.2 


4.516 


7.2+ 6.6 


6.381 


5.4+2.4 


8.904 


7.4+4.8 


16.570 


3.2+2.9 


16.518 


4.6+4.7 


1999 


2.030 


14.3±18.3 


5.133 


6.2+ 3.8 


5.649 


2.1+3.0 


8.776 


3.3+1.1 


16.851 


2.2+2.0 


16.470 


2.1+1.0 


2000 


1.647 


5.3± 4.3 


3.460 


4.3+ 3.5 


3.561 


4.4+4.4 


5.570 


4.2+3.5 


9.204 


4.2+3.8 


7.770 


3.8+3.6 


2001 


2.251 


5.5± 5.3 


2.906 


7.3+ 6.2 


3.671 


4.7+6.0 


5.220 


5.9+6.1 


7.249 


7.3+6.0 


7.760 


8.1+7.0 


2002 


3.686 


3.5± 0.4 


3.016 


3.0+ 1.3 


5.041 


2.7+2.4 


6.443 


4.9+5.3 


7.792 


3.6+1.4 


10.742 


4.1+5.1 


2003 


3.004 


5.7± 5.0 


4.186 


2.9+ 1.7 


4.447 


4.9+2.4 


7.560 


3.4+1.3 


9.104 


4.0+2.3 


14.034 


2.9+2.1 


2004 


3.261 


3.2± 0.8 


2.722 


3.3+ 1.8 


4.340 


6.1+2.3 


5.299 


1.9+1.3 


6.619 


5.7+2.4 


9.922 


3.9+2.1 


2005 


2.908 


4.3± 1.3 


2.350 


7.1+ 6.7 


4.122 


4.4+3.4 


4.687 


5.9+6.2 


5.756 


4.5+3.1 


8.575 


5.5+4.3 


2006 


2.906 


1.7± 0.4 


1.938 


1.7+ 1.7 


3.979 


3.0+2.7 


4.115 


3.0+1.9 


5.784 


2.7+2.8 


10.459 


1.9+1.3 


mean 




6.9±4.1 




6.0+3.4 




4.9+2.2 




5.6+3.0 




4.8+2.1 




5.2+3.0 



Year 




Hb(-5) 


Hb(+5) 




Hb(-4) 


Hb(+4) 




Hb(-3) 


Hb(+3) 




Flux 


(e+o-) 


Flux 


(e+o-) 


Flux 


(e+(r) 


Flux 


(e+o-) 


Flux 


(e+or) 


Flux 


(e+o-) 


1996 


2.511 


3.2+0.17 


2.535 


6.5+ 5.3 


4.965 


7.0+5 


4.470 


7.1 + 1.5 


8.003 


7.6+6.1 


5.838 


6.6+2.5 


1997 


1.303 


17.8+ 2.1 


1.507 


11.6+ 2.5 


3.134 


4.8+ 5.9 


3.813 


4.7+2.0 


6.391 


3.5+2.0 


5.381 


2.4+0.6 


1998 


0.923 


13.0+ 9.7 


1.735 


5.4+ 2.4 


2.393 


4.7+ 3.9 


3.318 


4.1+2.5 


5.268 


3.0+1.6 


5.648 


2.1+1.9 


1999 


1.215 


19.3+14.0 


1.594 


12.1 + 10.2 


1.674 


11.1 + 10.4 


2.786 


6.1+3.2 


4.403 


3.1 + 1.9 


4.322 


5.5+3.8 


2000 


0.255 


6.5+ 7.1 


0.822 


14.9+ 5.7 


0.626 


7.5+ 3.2 


1.593 


3.0+1.3 


1.654 


3.2+1.9 


1.977 


2.2+1.6 


2001 


0.642 


10.5+ 8.9 


1.068 


15.2+11.7 


1.064 


12.8+11.4 


1.620 


9.2+9.6 


1.923 


3.2+4.2 


2.210 


7.8+8.1 


2002 


1.253 


10.3+ 6.0 


1.011 


19.2+10.5 


1.684 


8.3+6.1 


2.217 


7.4+4.8 


2.165 


5.7+5.6 


3.235 


5.7+5.3 


2003 


1.000 


9.7+ 6.1 


1.568 


7.8+ 2.6 


1.515 


4.6+ 2.1 


2.715 


5.2+3.4 


2.627 


4.1+2.5 


4.347 


4.5+3.3 


2004 


1.161 


9.5+ 6.5 


1.427 


12.4+ 7.0 


1.568 


8.1+ 6.2 


1.861 


5.1+3.7 


2.062 


6.7+5.7 


2.841 


3.3+1.9 


2005 


1.062 


3.8+ 3.8 


0.969 


10.5+13.0 


1.458 


3.5+2.3 


1.535 


3.4+2.3 


1.818 


3.5+3.2 


2.609 


8.2+2.8 


2006 


1.012 


6.6+ 5.5 


0.886 


9.2+11.5 


1.536 


4.2+ 3.2 


1.323 


5.2+4.1 


2.024 


3.1+3.1 


2.841 


2.6+2.1 


mean 




10.0+ 5.2 




11.3+ 4.1 




7.0+3.0 




5.5+1.8 




5.2+3.0 




4.6+2.3 



Year 


Ha(-2) 


Ha(+2) 


Ha(-l) 


Ha(+1) 




Ha(0) 




Flux 


(e+cr) 


Flux 


(e+o-) 


Flux 


(e+o-) 


Flux 


(e+o-) 


Flux 


(e+o-) 


1996 


36.664 


9.4+6.8 


29.766 


11.0+7.0 


45.469 


9.1+6.8 


42.018 


11.6+7.6 


48.583 


11.1+7.6 


1997 


33.714 


7.0+8.2 


28.349 


8.1+7.1 


43.930 


6.8+7.1 


41.554 


7.9+9.2 


48.475 


7.4+8.0 


1998 


27.717 


3.5+3.2 


26.397 


4.6+4.0 


34.577 


3.5+3.8 


36.438 


6.4+5.5 


39.903 


7.3+4.9 


1999 


29.420 


1.7+2.1 


27.907 


1.4+0.8 


37.880 


1.5+0.1 


42.658 


3.5+2.6 


48.223 


2.5+2.2 


2000 


15.153 


3.9+4.2 


12.928 


4.4+3.7 


19.099 


4.6+3.9 


20.783 


6.6+5.7 


22.577 


7.5+6.8 


2001 


11.216 


8.1+6.3 


9.774 


8.0+6.2 


13.240 


10.4+8.4 


11.934 


10.1+8.4 


13.292 


16.2+13.0 


2002 


12.285 


4.1+1.1 


12.612 


6.0+2.1 


18.827 


4.7+3.8 


16.231 


14.9+8.3 


19.611 


12.4+4.7 


2003 


16.800 


2.5+2.2 


16.071 


4.6+4.5 


24.724 


3.6+2.5 


21.440 


13.5+9.0 


25.761 


9.8+4.9 


2004 


9.552 


4.8+1.5 


12.490 


9.1+3.2 


14.567 


1.4+0.4 


13.621 


18.0+4.9 


14.724 


9.7+2.5 


2005 


7.416 


4.3+3.2 


10.705 


9.8+5.3 


9.008 


8.3+6.2 


9.852 


32.0+21.3 


8.791 


26.0+25.4 


2006 


8.489 


2.6+2.9 


13.613 


2.6+2.0 


12.402 


4.2+1.6 


12.362 


9.3+4.9 


12.101 


12.8+5.6 


mean 




4.7+2.4 




6.3+3.1 




5.3+3.0 




12.2+7.8 




11.2+6.1 



Year Hb(-2) Hb(+2) Hb(-l) Hb(+1) Ha(0) 





Flux 


(e+o-) 


Flux 


(e+cr) 


Flux 


(e+o-) 


Flux 


(e+cr) 


Flux 


(e+o-) 


1996 


10.383 


7.0+4.5 


8.673 


6.2+1.5 


11.245 


6.3+4.5 


11.531 


6.7+2.2 


12.884 


7.8+3.9 


1997 


9.883 


3.7+2.3 


8.188 


3.3+3.1 


11.509 


4.1+2.5 


10.986 


4.0+3.1 


13.094 


4.9+1.1 


1998 


7.688 


2.7+2.5 


7.905 


2.1+1.4 


8.515 


2.7+2.9 


9.689 


1.7+1.4 


10.530 


1.8+2.1 


1999 


6.392 


4.1 + 1.9 


6.609 


5.1+2.9 


6.829 


6.3+3.9 


8.664 


4.8+3.2 


9.542 


6.6+6.8 


2000 


2.499 


2.6+2.0 


2.698 


1.7+1.0 


3.105 


3.9+2.5 


3.844 


3.0+1.7 


4.193 


5.5+4.6 


2001 


2.613 


4.4+2.5 


2.716 


5.6+5.6 


3.220 


4.5+4.0 


3.367 


2.5+1.3 


3.772 


8.0+8.5 


2002 


2.848 


5.0+4.2 


3.575 


4.9+3.8 


4.058 


5.3+4.2 


4.222 


5.3+4.2 


4.873 


6.8+4.6 


2003 


4.171 


4.2+2.1 


4.950 


3.0+1.5 


6.619 


2.9+2.9 


6.203 


3.4+3.0 


7.277 


4.2+3.2 


2004 


2.280 


5.4+2.5 


3.476 


3.3+2.7 


3.402 


4.6+4.6 


3.801 


3.5+3.1 


4.153 


3.9+5.3 


2005 


2.111 


4.0+3.4 


3.379 


5.0+2.4 


2.790 


4.4+2.5 


3.395 


7.7+2.7 


3.193 


9.1+4.1 


2006 


2.604 


2.5+2.5 


3.930 


1.2+0.9 


3.007 


2.0+1.7 


4.105 


2.3+1.8 


3.686 


3.5+4.0 


mean 




4.1 + 1.3 




3.8+1.7 




4.3+1.4 




4.1 + 1.9 




5.6+2.2 
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Fig. 13. The flux of different Hor segments as a function of the flux of the central Hor segment. The line flux is given in erg cm ^s ^ . 




I 



-5000 -4000 -3000 -2000 -1000 1000 2000 3000 4000 5000 

Velocity kms"' 



0.8 



0.6 I 



B 



0.4 



0.2 



Fig. 15. The variations of the Balmer decrement for different parts of the broad emission lines. The flux is given in erg cm ^ s ' A ' 



